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Abstract 

Effective field theory of interacting BFKL pomerons is investigated and Langevin equation for the 
theory, which arises after the introduction of additional auxiliary field, is obtained. The Langevin 
equations are considered for the case of interacting BFKL pomerons with both splitting and merging 
vertexes and for the interaction which includes additional "toy" four pomeron interaction vertex. 
In the latest case an analogy with the Regge field theory in zero dimensions (RFT-0) was used in 
order to obtain this "toy" vertex, which coincided with the four point function of two-dimensional 
conformal field theory obtained in |31j . The comparison between the Langevin equations obtained in 
the frameworks of dipole and RFT approaches is performed, the interpretation of results is given and 
possible application of obtained equations is discussed. 
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> 1 Introduction 

^ The scattering process of hadrons and nuclei in QCD with a large number of colors Nc have been 
■ - - vigorously investigated in recent years along two main lines. The Color Glass Condensate (CGC) approach 
[5l [6] , is formulated in the transverse position space and in large Nc limit the evolution dynamics of the 
model is analyzed in the terms of color dipoles [7]. Another approach, the QCD Reggeon Field Theory 
(QCD-RFT), was written and investigated in the momentum space for the BFKL pomerons, [21 [3l H], 
considering the picture of interaction of the pomerons in t-channel. The history of developing of such a 
picture begins from the pioneric paper [1] and it bases on the standard diagrammatic calculus developed 
for the interacting BFKL pomerons, having , as well, transverse position space formulation in large Nc 
limit of the theory, see [HI [9l UHl dH (El [13]. Both approaches are expected to be valid also in the case of 
symmetrical treatment of target and projectile and may include not only a vertex of pomeron splitting 
but also a vertex of merging of two pomerons in one. It gives the possibility for attempts to formulate 
and calculate the pomeron loops contribution to the scattering amplitude and makes these approaches 
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formally similar. But, in spite of this formal similarity, these approaches are written and formulated in 
the different frameworks and , therefore, the identity of the approaches is not fully clear. 

The QCD-RFT formulation of the high energy scattering uses a Lagrangian language describing the 
processes of scattering at high energy, [TUl ttH [121 13] , treating the processes of the high energy scattering 
in the terms of nonlocal effective field theory , which is based on the conformal invariant propagator of 
the BFKL pomerons and identical vertexes of pomeron splitting and merging. This approach has a long 
history of developing in the form of fenomenological RFT, [14^ I15j . and keeps a lot of mutual features 
with this "older brother", see [16]. The dipole CGC and so-called JIMLWK approaches are based on 
the consideration of the color dipoles as a main degrees of freedom in the high energy scattering in limit 
of large A'^c , O [7]. It is proven , that on the level of the zero transverse dimensions both approaches 
describes the same physics, [IT], in spite of the fact that these two approaches use very different pictures 
for the description of the scattering process. In QCD-RFT the calculations are based on the picture of 
t-channel propagating and interacting pomerons, whereas the CGC approach uses the picture of evolution 
of the dipoles of target and projectile in s-channel with interactions of dipole showers after the evolution. 
Therefore, considering the problem of pomeron's loops contribution into the scattering amplitude, both 
theories support very different strategies for the accounting of the loops. In QCD-RFT, because of 
the Lagrangian formulation of the theory, a program of such calculation may be formulated as usual 
perturbative calculations and these calculations in general must be very difficult, as it is always happens 
for the loops calculations in a field theory. In CGC approach the problem of the calculations of the 
pomeron's loops usually formulated as a problem of evolution equation with both , merging and splitting 
vertexes included and may be explicitly formulated in the terms of some effective Hamiltonian [18j. A 
Hamiltonian formulation of the theory has been considered in QCD-RFT as well, |13j . but the effective 
Hamiltonian formulation of the high energy scattering processes is out of the scope of the paper. 

In this paper we will investigate another possible way to account the pomeron loops in the different 
high energy scattering approaches, namely a approach when loops contribution are compactly coded in the 
Langevin equation. On the level of CGC approach the different types of arising Langevin equations were 
considered in [W\ . The correct impact parameter treatment of the pomeron loops was obtained 

in [20', '2T] , whereas the Langevin equation of [19] may be considered only as a some effective model of a 
correct Langevin equation. Indeed, in [19J the semiclassical treatment of large impact parameter behavior 
of the amplitude was used, when the sizes of the interacting dipoles are neglected in comparison to the 
impact parameter of the problem, that excludes a correct treatment of the pomeron's loops. Therefore, 
talking about the comparison of the Langevin equations obtained in the framework of QCD-RFT approach 
with the Langevin equations of the dipole approach, we will have in mind mostly the comparison of the 
results of this paper with the results of [20]. Nevertheless, we also will give the derivation of results of [19] 
in the framework of QCD-RFT . It must be noticed here, that the equivalence of QCD-RFT and CGC 
approaches on the semiclassical level was considered as well in the framework of a generating functional 
approach and BFKL pomeron calculus of \22\ [23] , where the problem of mutual description of QCD-RFT 
and CGC was formulated and resolved on the semiclassical level. But, as it was mentioned above, in our 
paper, in spite of the papers [El [23], the Langevin equation will be obtained with proper consideration 
of an whole impact parameter structure of the theory, where as a starting point the QCD-RFT approach 
will be used. 

The paper is organized as follows. In the next section we will introduce and consider the main 
construction blocks of the QCD-RFT approach and clarify they relations with such physically relevant 
quantity as unintegrated gluon density. In the Section 3 we will obtain the Langevin equation for the 
theory with only triple pomeron vertexes of splitting and merging. In the Section 4 we will introduce a 
"toy" four pomeron interaction vertex and will consider the possibility to write the Langevin equation for 
the theory with such additional vertex. The Section 5 will contain the comparison between the Langevin 
equations obtained in s-chanel dipole model and t-channel QCD RFT model. The Section 6 is a discussion 
of obtained results and conclusion of the paper. 
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2 Effective field theory of interacting pomerons 



In this section we will formulate the main results concerning the pomeron effective theory considered 
in [10\ nn [T2I I13j . A first ingredient, which we will need in the further calculations, is a action of the 
pomeron effective theory: 

S = So + Si, (1) 
where Sq and Sj are the free and interacting parts of the action correspondingly : 

So = /..d, .v,.^..v;dVi#fo,n.r.)G:i,,(n,..ir>;)*fe,.;,.;), (2) 



y-y 



and 



Si = / dy / j ^ % (^13 ^(y, n, ^3)) ^^(y, n, rs) ^2, rg) + (3) 



2as^c r , f d^rid^r2d^r3 . ,. ^, . ^, , 

H / dy / 2 — 2 — 2 — {Ln^^{y,ri,rs))^[y,ri,r2)^{y,r2,r3) . 

J J ^2^23^31 



In comparison with the action of [12j we omitted here the part of the action responsible for the interacting 
of pomerons with the target and projectile, the source terms, because for our following calculations these 
terms are not important. In the following derivations we also will not especially underline the fact that 
such quantities as and/or ki always denote the two dimensional vectors, it will be denoted as vectors 
only in the cases when precise definition of the vector structure of rj and/or ki will be needed. The part 
of the action, given by expression Eq. ([3]), reproduces the triple pomeron vertex in the large Nc limit 
with the use of an operator L13: 

Li3 = rj.plpl = rj.Vjvl . (4) 
The propagator of the theory, G^^^,{ri,r2\r[,r'2) , is defined throw the BFKL Hamiltonian [^[15] : 

G-\>{ri,r2\r[,r'2) = (^V^V? (^1- + H{rur2)^^ S^n - r[) 5\r2 - r'2)5{y-y') . (5) 

Now, let us consider only the triple pomeron interaction terms of the action and let us make a well known 
change of the variables in the action 



^(ri,r2) . '^Hn,r2) 
{ri,r2) 2 ,r{ri,r2) 2 ' 



(6) 

12 '12 



obtaining 



/■ , . f ,2„ ,2„ ,2^ ,^13 



Si = — ^ — - dy d^rid'^r2d^r3 {rl^(p{y,ri,r3))) (p\y,ri,r2) (j)\y,r2,r3) + (7) 



2a^JVc f f d^rid'^r2d'^r3 {r1i(f)\y,ri,r3,)))(p{y,ri,r2)(l){y,r2,r3) . 



IT 

In order to relate the amplitude with the unintegrated gluon density , as a next step, we perform Fourier 
transform: 

Hy,n,r2) = I ^^0^e-"^^^-"-^^^-'^^H{yM,qi-h) (8) 

for the functions 4>{y, ri,rj) in Eq. In the terms of the functions (p{y, ki, qi — ki) the triple pomeron 
term in the action will have the form: 

Si = ^"^^^ J dy J d'^k3d^q3d^k2{L3(t){y,k3,q3-k3))ci}^y,q3-k3,-k2)<PHy,k2,k3) + (9) 
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_^ 2a^_/Vc, I I (fk3(fq3(fk2 {Ls(t)^y,k3,q3 - k^)) (j){y,q3 - k3, -k2) (j){y,k2,k3) 



vr 

where, 

L3 = '^Ikl{q3-k3fvl, (10) 

see Appendix A for the detailed derivation, and as usual, q and k here are two dimensional vectors. Now 
we consider the case of interaction with zero momentum transfer, when: 

q3Hy7k3,q3 - h) = J q3Viy,k3)5'^{q3) = ip{y,k3). (11) 

In this case in Eq. ([9]) we have two delta functions arisen from the substitutions: 

Hy^qs - k3,-k2) ip{y,q3-k3)5^{q3 - k3 - k2) (12) 

kyMM - ^{y,k3)5^{k3 + k2) (13) 

that gives after the integration over k2 and (73 the triple pomeron vertex which is local in momentum 
space: 

Si = j dyj d\ (Vlk^ Vl^{y,k)) ^\y,k) ^{y,k)^ + (14) 

+ ^4^/ dyj d\{Vlk^ Vl^Hy,k))^{y,k)^{y,k). 

It is clear, that the function ip{y, k) is a scattering amplitude for the case when we omit a momentum 
transfer, taking it equals zero, i.e. for the case of forward scattering. In the impact parameter represen- 
tation this approximation corresponds to the semiclassical approximation of the large impact parameter 
limit, where the sizes of the interacting dipoles are neglected in comparison to the whole impact param- 
eter of the problem. This (p{y, k) function we may connect with the unintegrated gluon (parton) density 
function f{k) with the help of the following expression |1H I26j: 

f{y,k) = ^k''vl^{y,k), (15) 

and, therefore, the physical meaning of the (/^(y, k) became to be clear: through Eq. ()15p 99(2/, k) defines 
the unintegrated gluon density function for the processes of forward scattering. Basing on the Eq. Q- 
Eq. (llOp and Eq. ()14p -Eq. ()15p it is easy to generalize the expression Eq. (llSp for the case of non-zero 
momentum transfer: 

/(y, k,q-k) = ^k^ {q - kfvl 4>{y, k,q-k), (16) 

here (/)(y, k,q — k) is the amplitude defined by Eq. ([8]) and f{y, k, q — k) is a generalized (skewed) gluon 
(parton) distribution function. From Eq. and Eq. © it is easy to see, that our initial amplitude 
^ifi,T2) is simply Fourier transform of the f{y, k,q — k) and vice versa: 

f{y, k,q-k) Nc f d"^ n d^ r2 



/fe2 (q - ky 27r2 



/ (2 J' e-^^+-^(''-^)^(n,r,) (17) 



that determines very clear and transparence meaning for the amplitude $(ri,r2). Considering the equa- 
tions Eq. (jlSp and Eq. (jl6p . which are the Poisson's type equations, with the help of Green's function 
for the two dimensional Poisson's equation, we can write the inverse relation between f{y,k,q — k) and 
(f>{y,k,q- k) : 

4>{y,k,q - k,ei) = J d^k' g{k,k\q) ^^^^^^^7^)^ ^ 
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Figure 1: The relative positions of the angles for the k, k' and q vectors. 



^ / log p^ - 2kkCos{9i-e2) + k^ \ r f{y,k ,^2) 



J 47T ^\ A;2 J J {k'Y{q-k'Y ' 

The relative positions of the angles 9i and 62 are denoted in the Fig. [TJ Now it is easy to obtain the 
inverse expression for the Eq. ()15|) : 



^.(s,*:) = J d'gcl,(i,,k,q-k,ei=0)S'{g) = (19) 



J ^ k"^{q-k'Y J 47r k"^ ) 



A^c 7o 47r Jk^ k'^ Air 

see also ^6]. With the help of Eq. (|14p and Eq. (jlOp we could obtain the expression for the triple pomeron 
vertex in terms of unintegrated gluon density, see [26l|l6] and Appendix B. In general, using Eq. ([7]) and 
Eq. (jlSp . it must be also possible to obtain the Bartels triple pomeron vertex written in the terms of the 
skewed unintegrated gluon densities f{y, k,q — k) , see but we do not consider this task in this paper. 
So, in the further consideration, using the amplitudes <^(y, rj,rj) or (j)(y,ri,rj) defined in the transverse 
position space we will always remember, that these quantities are related to the generalized (skewed) 
gluon density function f(y, k,q — k) in the momentum space. 

3 Langevin equation in the theory with the triple pomeron vertex 

3.1 Langevin equation for the $(?/,rj,rj) field in the transverse position space 

In order to introduce an auxiliary field in the theory we come back to the particular part of the action 
from the Eq. ([3]): 

= f dy f f (Li3«>(y,n,r3))cl>t(y,ri,r2)c^t(y,^2,,3) (20) 

and we rewrite this expression in the following form: 



•5$ = / dy 4 ^\y,Pi^P2) / dy 3 $'(y,/92,/0i 



(21) 



n4 ^4 



Now we introduce on the scene an auxiliary field through the Gaussian integration over the auxiliary 
field V' : 



^ F12 \ ^12"l'2'^ll' 



e^'' = N I D[i;]exp{-- I dy \ i^(.y,Pi,P2) — {L^^' Hy, pi, p^)) 222 



-V)5'^{P2- P2) [ c^y' ^-T— ^ V'(y',/'2>Pi) - / dy il;{y,pi,p2) <^\y,pi,p2) \ , (22) 

where a distribution functional for -0 has the following form: 



^1 



^ P12 \ ^ P^2Pi'2 Pll' ) 

5{y-y)5\p2-P2) I dy'^^^^i^iy\p[,P2)] , (23) 

J /'l'2' J 

with ^ 

AT = (^1 . (24) 

Let us define a new part of the action of the theory: 

^ If, d'^pi(fp2 ^fo^l^c,^ ^, /'l2Pt'2' /o^.^ 

= -77 / dy 5 ij{y,pi,p2) — (L^-^f ^{y,pi,p^)) — — ^ — (25) 

-I Pl2 \ ^'^ '^12^l'2^1l' / 



Ky -y)^\P2- P2) I dy ^ Tpiy ,p2,Pi) 



I ,0 / 

/ / / . 



^1'2' 



and rewrite whole action in the following form: 

5 = So + Saux + 5'$t - / dy ^^-^^^-^^ ip{y, pi, P2) ^\y, pi, P2) , (26) 

J P12 

where as a S^'f we denoted a second part of the action Sj from Eq. ([3|) . Writing the equation of motion 
for the ^{y,pi,pj) field 

S S 

s^Ky,Pi,P2) = ° ^^^^ 

we see, that we obtain the equation for the field <I> which does not depend on the field <I>^ but rather on 
the field ip: 

G-'Hy,P„Ps) + ^ / -^^^Hy,PuP2)Hy,P2,Ps) (L13) - ^^^^ = 0, (28) 

^ J P12P23P3I Pl3 



or 



+ H{pi,p3)] <^>{y,puP3) + (29) 



2a^A^c /■ d'^p2p\^ , , , N / X 1 / 

^ / -^^^{y,PuP2)Hy:P2,P3) - {Li3)-'i^{y,Pi,P3) = 0. 
^ P12 P23 
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The field ^p in Eq. (j29p may be considered as a noise field with the following properties: 

< V(y,Pi,P2) > = ; (30) 



/ / / , 



< '^{y,Pi,P2) ,i^{y ,P2,Pi) > = ^"'^ ^^^^^ ^ {Lii' <^>{y,pi,p'i))6{y - y')5^{p2 - p'2) (31) 



that follows from the form of the distribution fmictional Eq. (j23p for the ■0- The last term in Eq. ()29p may 
be rewritten with the help of the following properties of initial Green's function, i.e. Green's function at 
zero rapidity. Indeed, let us consider this Green's function Go{pi, p^lp^' , p^'), which has a form, see [25] . 

22 22 

Go(pi,p3|Pi',P3') = In In (32) 

and satisfies 

(pi3)~^Li3 (Go(pi,P3|/3i',P3')) = (2vr)'^(52(pi -/?!') (52(P3 -/>3')- (33) 
Using Eq. ([33|) we obtain an operator identity: 

(V^J-^ (V^3)-^ = / ^^^^Go{p^,p^W,p,>). (34) 
With the help of Eq. (I34p we rewrite Eq. ()29p in the following form: 



(^-^ + H{pi,ps)^ ^{y,pi,P3) + 
2alNc f (fp24i , . , . r (Pp^jfp^ ilj{y,pi,p3) 

/ — —'^{y,Pl,P2)'^{y,p2,P3) - / x4 Go{pl,P3\pi',p3') 4 = 0. 

vr J Pto Pr>^, J (^T^j Pis 



(35) 



P12P23 J V^^J FIS 

with the correlator for the function Tp{y,pi,p3) given by Eq. (fSTl) . Both Eq. (f29l) and Eq. (f35]) with the 
auto correlator for the noise field Eq. (j3ip are the required Langevin equations of the QCD-RFT approach, 
valid for any value of the impact parameter of the problem. These equations are pretty complicated, 
therefore, let's try to simplify them using the conformal basis representation for the pomeron fields and 
assuming special properties of the theory at high energy limit. 

The following property of the field ^{y, pi, pj) at high energy limit of the theory could help us to 
simplify the Langevin equation. Let us expand the ^{y, pi, P2) field on the conformal basis formed by 
the functions £^^(„,,.) ,po(pi, P2) , iZJ: 

E,in,.,,M^P2) = (-^) (/|r) (36) 

\PlOP20/ \Pi0P20y 

and 

Ll3 E^(n,u) ,poiPl^P2) = \,ln,v) ^t^M ,poiPl^ P2) , (37) 

with ^ 

K{n,u) = Ap = ((^+1)2 +4^2) ((^_1)2 +4^,2) ' (38) 

where we used the same notations as in [12 1. This expansion has the form 

Hy ^ Pli P2) — Efj,(n,u),po {pi,P2)<^p{y) = E E,,{pi,p2)^p{y) (39) 
fj. p 

where 



^/.(y) = / ^\ ^' E;ip,,p2)^y,puP2), (40) 
P12 
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see |27) . In Eq. ()39p and in the following expressions the notation of the conformal summation always 
means 

OO 2 TL^ 

E = E jdu^-^jd-po. (41) 

Ai n=-oo /I J 

At high energy limit , as it obtains for the BFKL equation at high energy limit, see for example [29], we 
can assume that the main contribution in the sum in Eq. (|39p comes from the minimal conformal weight, 
namely when n = and u = {). In this case we have 

Li2$(y,Pi,P2) = Li2 E^{pi,P2)^^{y) = E V^^/^(/'i'/'2)^/.(2/) - (42) 

~ E ^M(pi.P2)^M(y) = Hy,Pl,P2)- 

Here we used the fact, that in the high energy limit when n = and = we have X~^^^q ^^q^ = 1. The 
same approximation was used, for example, in [23]. It must be clear, that this high energy approximation, 
which is based on the behavior of the BFKL amplitude, i.e. single pomeron, may be not correct in general 
in the theory of interacting pomerons. Now, supposing that the Eq. (j42p is valid at high energy limit, we 
rewrite the Eq. ()35p and Eq. ()3ip in the following form: 

d \ 

— + H{pi,p3)j $(y, pi,P3) + (43) 
2a^^Nc f d^p2pli^, , f(fp-.>(fp^' 

/ — — —^[y,pi,P2)'^{y,p2,p3) - / — /n N4 Go[pi, pslPi' , Ps') i^{y, Pi, P3) = o. 

^ -J P12P23 V^^J 

< ilj{y,pi,p2) > = ; (44) 

< ij{y,Pi,P2) ,ipiy' ,p2,p'i) > = ^ <i>{y,pi,p'i)6{y-y')6'^{p2- p'2), (45) 

^[Pi2pI'2' Pn' 



were we made the following substitution tp{y, pi, P2) i^iy-, Pi-, P2) / Pi2- The equations Eq. ([^S^ -Eq. ([^5]) 
may be considered as a high energy approximation to the correct Langevin equation Eq. (|29p with the 
auto correlator Eq. (f3T]l . 



3.2 Langevin equation for the Lp{y,k) field in the momentum space 

Let us now consider the part of the action given by Eq. (I14p for the case of the forward scattering and 
formulated in the momentum space: 

=^4^ I dyj d-'k (Vlk^ Vl^{y,k))^Hy,k)^iy,k)^ . (46) 

As it was mentioned previously, a dual expression of the S^p part of the action in the coordinate space will 
represent a semiclassical approximation of the initial action Eq. ([T])-Eq. ([2]), and , therefore, the results 
of this subsection are related to the Langevin equation of ITU]. So, performing the same transformations 
as in the previous section, we rewrite the part of the action in the following form: 

S'P = lJdyJ d^k^\y,k)(^^{Vlk'vl^{y,k))6{y-y')6\k-k')\ J dy j d^ k' ^\y' , k ) . (47) 
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Introducing an auxiliary field ipiu^ k) we write for the c^^: 



6iy-y')6\k-k') J dy' cfk i^iy' ,k') - J dy d^ki^iy,k) ^\y,k) j . (48) 
Reabsorbing the (^ ^"^^^ (V| k'^'\l\ip{y, k)) ^ ^ factor in the definition of the field -0) we obtain: 
e^f = N J L»[V']exp|-^ J dyd^k'4){y,k) J dy' d^k ip{y' ,k') 5{y - y')5^{k - k') 



dyd^kij{y,k)ipHy,k)^^^iVlk^Vlip{y,k)) } . (49) 
As in the previous case, N here is 

N=(^J D[i,]expl^-^ J dyd^ki^{y,k) J dy' d^k i;{y' ,k) 5{y - y')5^{k - k')^^ ' . (50) 
The whole action for the ip{y, k) and ^p^y, k) fields now takes the form: 



S = So + Sau. + - J dyd^ki;iy,k)^Hy,k)J^^^{Vlk^Vl^{y,k)). (51) 



where 



Saux = -\ I dy I d^ki>\y,k). (52) 



Writing an equation of motion for the field ip{y, k) 

6S 



6(p^y,k) 

we obtain equation similar to Eq. ([28]) : 



^ — + H{k)j ip{y,k) + -^^{y,k) - (Lfc) U{y,k)^^^{Lk^{y,k)) \ = 0, (56) 



where correlators for the auxiliary filed ip{y, k) have the form: 



(53) 



Lk[-^ + H{k) ^{y, k) + (^(y, k) - Hv, k) J ^^(Lfc y^{y, k)) = , (54) 

\oy / TT V ^ 

with the operator Lk from the Eq. (1140 : 

U = Vlk^Vl. (55) 

We rewrite Eq. ([5lll : 



< ^iy,k) > = ; (57) 
< i;{y,k),4^{yi,k,) > = 5{y - yi) 5\k - ki) . (58) 
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It is important to underline, that contrary to the simphfications obtained for the equations Eq. (j43|) and 
Eq. (I45p for the ^{y, pi, pj) field at high energy limit, here we cannot write that Lkip{y,k) = ip{y,k). 
Indeed, such condition leads to the very non physical restriction on the unintegrated parton density 
function f{y,k) which follows from the Eq. (llSp and Eq. (llOh : 

Vl/(.,fc) = ^/J<ifc'^^log(^). (59) 

Of course, such identity can not be satisfied in general for arbitrary /(y, k) function. More of that, 
due the use of semiclassical approximation in derivation of this result, it is not clear at all, how such 
approximated Langevin equation may be used as a guideline for the pomeron loops calculations. Indeed, 
by definition, a loops calculation must correctly treat the impact parameter structure of the theory , 
whereas the Eq. (j56p and Eq. (j3.2p were obtained in the approximation of the large impact parameter. 



4 Langevin equation for the theory with "toy" four pomeron vertex 

4.1 "Toy" four pomeron vertex in the effective theory of the interacting pomerons 

The Lagrangian of the RFT-0 model for the q and p pomeron fields, which includes also a four pomeron 
vertex, we could define in the following form: 

L = qp + pqp - Xq{q + p)p + X q^ p^ , (60) 

see |141I15( [T7]. where is a bare pomeron intercept, A is a vertex of triple pomeron interactions and A a 
four pomeron interaction vertex. For the case of "fine tuning" of the vertexes, when t = tt, that defines 

a "magic" value of the four pomeron vertex A , the Hamiltonian of the problem has a factorized form in 
the terms of q and p fields : 

-H= p{q - -q^)p - X{q - \q^)p' = p{q - - q^) (p - -/). (61) 
p A p p 



We are not interesting in the further investigation of the RFT-0 model here, see more details in 117 ^ 123 ^ 130]. 
but as a guideline for the derivation of our "toy" four pomeron vertex, we will take the same as in Eq. (j6ip 
property of the factorizability of the Hamiltonian for the case of "magic" value of the four pomeron vertex. 

We will begin from the free part of the effective pomeron theory action, Eq. ([2]), written in the 
following form: 



+ 2 



^ ^ 1 y ^yl^^i^^^s /^^t,^_ ^ , 5(Li3$(y,ri,r3)) 5(Li3$t(y, n, rg)) 



'^y— \^\y,ri,n) — — $(y,ri,r3)l + (62) 



dy 4 ( ^>(y,ri,r3)(Li3i7(ri,r3) (y, n, r3)) + (Li3i?(ri, r3) $(y, n, r3))$^(y, ri, r3) 

'^13 

Here H{ri, r^) is a BFKL Hamiltonian, [ 2411251 [27 J. In order not to confuse this Hamiltonian with general 
Hamiltonian of the problem H, the BFKL Hamiltonian and operators related with BFKL Hamiltonian will 
be always written with the arguments of the BFKL Hamiltonian, namely in the form H{ri,rj). Writing 
the Hamiltonian of the problem in the conformal basis, formed by functions E^[ri^r^) of Eq. (|36p . we 
will omit in the further expressions the common integration factor 

ivl^^. (63) 

J ^13 
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The expression for the " free" part of the Hamiltonian in the conformal basis is the following: 

^^0 = \Y.E^{n,r^)u:,X-^'^,{y)Y,E:[rur^)^l{y) + (64) 
+ ^E^M(^i'^3)^M(y)E^^(^i'^3)^.A;i<I>t(y), (65) 

[1, V 

where we used Eq. ([37|) and 

H(rx,r-i)E^,{rx,r-i) = - ^^(n, rs) , (66) 

with the eigenvalues for the eigenfunctions Ep,{r\^r'i) of the BFKL Hamiltonian [28]. Before the 
definition of the "interacting" part of the Hamiltonian, let us introduce functions 



/dp' To 
^-^$(y,ri,r2)$(y,r2,r3) = ^ E^(ri, rg) ^'^(y) , (67) 
'^12 ^"23 



and 

^^(y,r-i,r3) = / ^^$^(y,ri,r2)$t(y^^2,r3) = E £;;^(ri,r3)^'t(y). (68) 

J '^12 ^^23 ii 

With the use of the ^(y,ri,r3) and ^^(y,ri,r3) functions the "interacting" part of the Hamiltonian, 
which corresponds to Eq. ([3]), obtains the following form: 

= ^3) A;1 ci>^(y) ^ ii;:(ri, rs) (y) + (69) 

+ ^^E^;(^i'^3)A^ic^t(y)^i?.(ri,r3)^.(y). (70) 
Let us now assume the following anzats for the action with the "toy" four pomeron interaction vertex : 
S,p = -C f dy f f^:ip:^F,p E,{n,rs)^,{y) J2 E:{n,r,) ^Uy)] = (71) 

= f dy f fllpl (Y^{F,E,{n,n))^,{y)Y,{F2El{n,n))¥M + 

^13 V ^ u 
+ ^ (F2 E^{n,T^)) ^,{y) J2 (A E:{n,r^)) ^l{y) ) = 

fj. V j 

= f dy f fllpl h^E,{n,r,)mMT. h'^K{rur,)^l{y) + 

+ E /2M^M(ri,r3)M/.(y) E /i. ^^(n, rs) *t(y) ) , 

M V J 

where we introduced some operators F\ and F^ such that 

F,^^(ri,r3) = /i^S^(ri,r3). (72) 
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Finally, again omitting the integration over J dy J (PridPr^ / rf^, we write the Hamiltonian which is 
corresponding to this anzats: 

H^P = CY. fi^E^iri,r^)^,{y) ^ /a. ^^(ri, rg) (y) + (73) 
+ CY.h^E^{r^,r^)^MY.hvEl(r^,T^)^l{y), (74) 

where the constant C will be defined in the further derivation of the Hamiltonian. Now, collecting all 
terms Eq. ([Ml) . Eq. ([69]) and Eq. ([73]) together, we obtain: 



fl u \ 



+ (75) 



^.(2/) - (76) 



- ^T.E,iri,rs)^,iy) ^ Kiri,r3)K' fe(y) - Xufi,f2. ^tiy)) - (77) 

- ^E^;(^i'^3)^';.(y)E^-(^i'^3)A;^(«>.(y)-^A./i^/2.M/.(2/)) . (78) 

fl u \ s c / 

Assuming for Eq. (|75|) the same factorization property as for the Hamiltonian of RFT-0, we find the 
following values for the C , /i^ , f2u'- 

C = 2 {^^^ , fl, = , f2u = ■ (79) 
Continuing the derivation, we obtain for the Hamiltonian: 

H = E^M(n,r3)E^:(-i'-3)(-.A;i^-4r^^.(y)A;i) (80) 
(^^tiy) - + (81) 

+ E^;^(n'-3)E^^(n,r3)(o.,A,^^-^*;.(y)A;^) (82) 

($.(y) - ^^;^^.(y)) . (83) 

For the further simplification of the expression Eq. (j8U|) we will do the following. First of all let us 
consider the , ^'J, functions in the Eq. ()80p . From the Eq. ()67p we have : 

(fr2 rf 

fi.,v-' '12 '23 fji 

and using the orthonormalization properties of -E^, [2^, we obtain: 

*M(y) = E / ^^^|^^^4^^«'(n,r2)i?.(r2,r3)i?;(ri,r3)cl>^(y)cl>,(y) = F^,^,. $^(y) <I>.(y) , (85) 

u),!/ ''12 ^^23 ''13 



E / ^^i?/.(n,r2)^.(r2,r3)cl>^(y)'3>.(y) = E E^{r^,T^)^ ^{y) , (84) 
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where V^^ui.^ is the triple pomeron vertex in the conformal basis, see [31j, and summation over repeating 
indexes w, v is assumed in the form of Eq. (I4ip . Another observation is concerning the omitted integration 

/ '^^ . This integration over conformal functions E*{ri, r^) £',y(ri, ra) in the Eq. ([50]) gives 5uu, that, 

'"13 ^ 
excepting the integration over rapidity, determines the full Hamiltonian: 

(86) 

in full analogy with the RFT-0 Hamiltonian Eq. (|61|) . We see now, that the "toy" four pomeron vertex 
in the conformal basis has the form 



v,P = (^4^) E ^M^'^{y)v-,,^,uV^^^',,'^l>{y)^l>{y) 



and the same vertex for the action in the usual field basis looks as follows: 



(87) 



- , , 111 1 

TT I J J J rf2r^3r^2'V3 



+ (-H'"^(ri,r3)^>(y,ri,r2) $(y,r2,r3)) Li3$^(y, n, rg) ^>^(y, , r3)) | . 

From the form of this vertex it is clear, that at least for this "toy" four pomeron vertex in the usual 
basis, the procedure described in the previous section gives very complicated auto correlator for the noise 
field. Indeed, now, due to the very complicated form of the kernel in the Eq. (|88p . this auto correlator 
will include a complicated expression with the L and operators with the square of the ^{y,ri,rj) 
field. Therefore, we will write this expression only in Appendix C. Nevertheless, for the action in the 
conformal basis the possible expression for the auto correlator of the noise field looks much simpler and 
we will consider this derivation in the next section. 



4.2 The Langevin equation for the "toy" four pomeron vertex in the conformal basis 

The action of the effective pomeron field theory in the conformal basis looks as follows: 

-1 anHy) 



dy 



-u;,X-^' (^^,{y) - ^^u;~%,^,My)'^Ay)^ (^^Uv) " ^^^^^ V^,^,,cl>t„(y)cl>t (y)^ | , 

where we used the kinematic part of the action from the equation Eq. (I62|) . The action describes the 
interaction of infinite number of the pomerons through infinite number of different triple pomeron and 
four pomeron vertexes in accordance with the results of [12^ [T3] . Let us now write the part of the action 
which contains the square of the field: 
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Proceeding as before we obtain for the 5$: 

.2 



\ ^ y ^ 

and for the auxihary field action we have: 



K'u:-^^^^,{y)<^^,{y)V-,,V,,^A 5{y-y') / dy ^l{y') (91) 



Saux = ~\ j dyil^w{y)\^ ^ 4a^iV^ \^^^^{y) V^^ui^u - 

2(^4^1 K'^~,'^..'^y)^Ay)y',,^',u'y,,^,^] My)- m 



vr 

Now, whole action Eq. (I89|) may be written as 



s = j dyY.U%^y)K —d^-2^^^y^^^ —t- 



-u:^\-^^^{y)^l{y) + ^£-^A^i$^(y)c&,(y) y^,^,,cl>t (y) | + Sau. - J2 J dyMy)'^^- (93) 
The equation of motion for each <I>^ field in this case has the form of Langevin equation: 

d^^iy) _ 



dy 



- ^M^u{y) Vf,,^,^ + , (94) 



where we are not summing up over the index /i. The field in Eq. (j94|) we can consider as a noise field 
with the following correlators: 

< tl;^{y) > = 0; (95) 

<i^M,Myi)> = ^^EV (96) 



<i>M(y) W-^4^^;''J'«,'(y)'J'.'(y)^;iyyW) (97) 



We see, that we obtained the infinite number of Langevin equations for the fields with very complicated 
non diagonal auto correlator for the noise field ip^. It must be underlined , that the evolution equations 
for the fields are also very complicated , there we have summation over w, v indexes in Eq. (j94p that 
means integration and infinite summation due the definition of the summation procedure in Eq. (I4ip . 



14 



5 The Langevin equation in the dipole approach 



Now we will consider the results on Langevin equation obtained in [20] in the framework of the dipole 
CGC approach and will reproduce our Langevin equation basing on the equations of [20j. Let us consider 
the Langevin equation of |2Qj for the process of the scattering of an arbitrary numbers of dipoles off 
the target. The noise (fluctuation) term of the Langevin equation, corresponding to this process, was 
obtained in [20j and it has the following form: 

f{ri,r2,y) = C I (fpi (fp2 d^P3 Go{ri,r2\pi,P3) Jvj^ Vj^<^{y,pi,p2)i^{pi,P2,P3,y) (98) 

J Pl3 

where we adapted the notations of [20] on the notations of the paper, here we introduced a <l>(rj, r^) field 
as a RFT counterpart of the Tyiri^Vj) function from [20]. The constant C in the expression Eq. (|98p is 
related to the possible difference in normalization of the amplitudes in both papers and we do not fix it 
in the expression Eq. pS]) . for our following derivation it is not important. The noise field i^lpi, p2, P3,y) 
in Eq. (IDSp is a noise field from [20] with the following auto correlator: 



/ / / / ^ 



< T^{pi,P2,p3,y) , i^{pi,P2,p3,y) >= ^'^{pi- P2)^^ip2- pi)s'^{p3- p3)^iy -y) ■ (99) 

Using Eq. ([99|) . we can define the auto correlator for the field f{ri,rj,y): 

<f{ri,r2,y),f{r[,r2,y)>=C^J d'pi d^p2 d^P3 ^ Go(ri, rsjpi, ps) V^V^, ^y, pi, pa) (100) 

f d^p[d^p'2d^p'^'^-^^Go{r[,r'2\p[,p'3) Jv^, V^, ^{y, p\, p'2) < u{pi, p2, p3,y) , u{p[, p2, p's,y') > . 

J P^l^l V ^2 
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The straightforward calculations gives with the help of Eq. (I99p : 



/ / / , 



< f{ri,r2,y) , f{ri,r2,y) >= C / d pi d /)2 d Ps G'o(ri, r2|pi, pa) G'o(ri, r2|/92, Ps) (101) 



P12 



^y, PUP2) JVl ^y, p2, pi) 5{y - y 

Pl3 P23 



Now we will use Eq. ()33p in order to rewrite Green's functions Go in the following form: 

Go(ri,r2!pi,p3) = V;^^V;^^6^{ri-pi)6\r2-p3) (102) 

and 

Go(r;,r;|p2,p3) = V-2v/52(r;-p2)52(r;-p3). (103) 

Liserting Eq. (jl02p and Eq. ()103p into the Eq. (jlOip it is easy to see, that operators and 
now may be reabsorbed in the definition of our initial field /(rj, r_,-, y): 

f{n, r„ y) ^ /(r„ r,-, y) . (104) 

With this new noise field f{ri,rj,y)^ the equation of motion for the <l>(rj,rj) field, or for the Tyiri^Tj) 
amplitude from PU], obtains precisely the same additional noise term as in Eq. (|43p . Indeed, now the 
auto correlator for the f{ri,rj,y) has the form: 

<f{n,r2,y),fir[,r'^,y')>=C^I d^p, d^p2 d^ps y^V^, V^^ $(y, pi, P2) .Jvj^ V% <I>(y, p2, pi) (105) 
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Pl3 P23 

Performing integration and interchanging variables r'l and r'2 in the /(r^,r2,y'), we obtain: 

2 

<f{n,r2,y),f{r'2,r[,y')>= 4^ (vl^Vl , ^y^n^^)) S^r^ - r'^) 5{y - y ) . (106) 
Comparing this expression with the auto correlator given by Eq. (|31|) : 

< i^[y,pi,P2) ,w{y ,P2,Pi) >= 2 — i^n' ^{y^pi^Pi))^vy -y )^ (p2 -P2)> (io7) 

we see, that if we reabsorb l//ofj coefficient in definition of the noise field Tp{y,pi,pj) in Eq. p07p . then 
we obtain the same as Eq. (I106p auto correlator for the noise field. 



6 Discussion of results 

As a first main result of the paper we consider a Langevin equations Eq. (pU]) -Eq. (f^T|) . Eq. ([5B]) -Eq. (j3.2|) 
and Eq. ()94p -Eq. (I95p obtained in the framework of QCD-RFT. The meaning of these equations is simple. 
The derivation of them is achieved with the use of the classical action of the theory, i.e. action without 
pomeron loops . The introduction of the auxiliary field in the classical action instead the square of the 
field leads to the Langevin equation for the field Iterating calculations procedure for the $ field 
together with the calculation of the auto correlator for the noise field tp, see for example discussion in 
|19j . will lead to the account of the pomeron loops contribution in the pomeron field Proceeding, we 
will calculate an all loops contribution to the pomeron field and will resolve whole quantum problem for 
this effective theory. Of course, practically, due the very complicated form of the auto correlators for the 
noise fields, this task looks at least not easy. 

The interpretation of the Langevin equations in this paper and the interpretation of the Langevin 
equations from the papers |19tl20j are the same, in spite to the different forms of obtained equations. As in 
|19j and [20j we also obtained two different Langevin equations formulated in the different frameworks. Let 
us , therefore, discuss the similarities and differences in the forms of the Langevin equations, formulated in 
the momentum and coordinate spaces separately. First of all, we consider the obtained Langevin equation 
Eq. (|56p for the (p{k, y) field in the momentum space. The Eq. (|56p looks very similar to the Langevin 
equations of [19]. There is only one move which reduces Eq. (I56p to the corresponding equations in [19] 
and this move is an assumption about the action of the operators L and of Eq. ([J]), on the function 
y). If we assume, that the (/^(A;, y) is a eigenfunction of L and with the unit eigenvalues then the 
corresponding equations will be the same. But, as it was argued in this paper, this move is impossible 
in the framework of the present approach because in this case a very unphysical restriction Eq. ()59p on 
the form of the unintegrated gluon density function /(y, k) is arising. Therefore, in general, without 
some approximation it is impossible to write Eq. ()56p in the form obtained in p^. More important, 
that the function (^(fe, y) is defined for the case of zero momentum transfer, i.e. there is no correct 
impact parameter dependence account of the amplitude is performed, and, therefore, it is in principal 
impossible describe loops contribution in the pomeron field with the use of the (^(/c, y) field. The Langevin 
equation for this function in the form of Eq. (|56p . therefore, may plays a role of some toy model which 
probably may describe some properties of real QCD. In general, this function also may be applied for the 
semiclassical solution of the problem, i.e. for the calculation of the "tree" pomeron structure, see |16j . 
Another application of this amplitude and semiclassical approach was developed in [22^ |23] , where the 
calculations of the amplitude were performed in the framework of a generating functional approach which 
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has a strong relations with the probabihstic interpretation of the BFKL pomeron interactions. It must 
be also noticed, that the approach developed in [23j has a mutual roots with the approach considered in 
this paper, focusing, nevertheless, mostly on the semiclassical solutions of the problem. 

Concerning the correct Langevin equation for the ^{y,ri,rj) field in the transverse position space 
we see, that the situation here is more obvious. The calculations in the previous section show the 
equivalence between the Langevin equation obtained in [20j in the framework of the s-channel dipole 
model and Langevin equation given by Eq. (j29p -Eq. (j3ip which was obtained in the framework of the 
t-channel QCD-RFT. In spite to the differently written forms of the noise terms and auto correlators of 
the noise terms , the Langevin equations in both approaches are the same. The Langevin equations it is a 
another description of the dynamics of the physical problem, and the same form of the Langevin equations 
in both approaches means a equivalent description of the quantum fields dynamics in QCD-RFT and 
dipole approaches. Using Langevin equation as a bridge we can pass from the QCD-RFT side to the 
dipole side of extended Balitsky-JIMWLK hierarchy and vice versa. We can conclude, therefore, that the 
ry(rj,rj) function considered in [20j in the framework of the extended Balitsky-JIMWLK hierarchy of 
equations is equivalent to the pomeron field considered in the framework of t-channel QCD-RFT, where 
in both cases all pomeron loops contribution is included in. More of that, due the equivalence between 
the QCD-RFT approach and dipole model of [20] prooven in this paper and due the fact that in [20] was 
established the equivalence between approaches of [20] and [21], we can conclude, that our QCD RFT 
t-channel model is equivalent to the modified version of the JIMWLK equation of |21) as well. On the 
level of RFT in zero dimensions this fact was clarified in ^7\, and now we see, that the s-channel dipole 
and t-channel RFT approaches are indeed equivalent in the physical space of two transverse dimensions. 
This result we consider as a second main result of the paper. 

The consideration of the theory with the "toy" four pomeron vertex in the conformal basis is similar, 
as it must be, to the derivations of |121ll3j. As in [12], we obtained the reformulated theory with infinite 
number of one dimensional pomerons, with only additional terms due the four pomeron vertex. In spite 
to the redefinition of the degrees of the freedom of the theory, the corresponding equations are not 
simpler then the equations in the usual basis. Indeed, the equation of motion of this theory contains 
infinite number of one dimension pomerons with an infinite number of different and very complicated 
interaction vertexes which include three integration and three infinite summation over repeating indexes. 
The equations in the form of Langevin dynamics also include very complicated correlators of Eq. ()95p . 
where the auto correlator for given conformal weights /u , v depends on all other conformal fields ^ of the 
theory. Nevertheless, as it was mentioned in [12j, we can truncate the sums on some value of conformal 
weight n and try to solve simplified truncated theory. Still, this task will be not easy and we leave it for the 
future studies. Concerning the form of the expression for this four pomeron vertex it is also interesting 
to note, that obtained in this paper vertex is the same as the four point function of two-dimensional 
conformal field theory derived in [31] from the conformal bootstrap point of view on the BFKL pomeron. 
In the given context it means, that conformal bootstrap in application to the interacting BFKL pomerons 
will lead to the factorized Hamiltonian of the problem with precisely zero ground state, see beginning 
of the Section 4. We do not know about the deep physical reasons for a such relation between the 
conformal field theory in two dimensions and the Hamiltonian of the interacting BFKL pomerons, and, 
therefore, we do not consider here this subject, which potentially may be very interesting. Finishing the 
theme of four pomeron interaction vertex, we can also notice, that our vertex is not the same as the four 
pomeron vertex introduced in the [22] on the base of the physical reasons and arguments different from 
the considered in this paper. 
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Appendix A: 



Let us consider the following part of the action: 



Si = I dy I ti^ri (i^r2 d^rs (r|i0"f(y, ri, rs))) n, r2) (/'(y, r2, rs) , (A.l) 



where we omitted the coefficient in front of Eq. ()18p . Let us make the Fourier transform of the (j) functions: 
Hy,ri,r2) = I ^^l^^e— ^(^^-^^)<^(y,A:i,gi-A;i), (A.2) 

</'(y,r2,rs) = I '^2fy' ^''-''^ky,k2,q2-k2), 

Operator ri^ throw the Fourier transform is changed to : 

''31 

¥ rl = rh V? r?s - K ^1 (^3 " ^s)^ ^l, = U ■ (A.3) 
Inserting Eq. (|A.2p and Eq. (jA.3p back in Eq. (|A.ip we obtain: 

Si = dy d^ ki d^ k2 d^ / d^qid^q2d'^q3{L3(f)'f{y,ks,q3-h))<P{y,k2,q2-k2)4>{y,ki,qi-h) (A.4) 



(i ri d r2 d rs ^^^^ (gi - fci) g-ir2 A;2 -ira (52 - ^2) gi^s ^3 (qs — ^3) 



(27r)2 (27r)2 (27r)2 
Let us consider the last line of Eq. (|A.4p : 



d ri d r2 d ^.^^ _j ^.^ _ ^j^) -ira (92 - ^2) girs fca +iri (93 — fca) (A 5) 



(2^)2 (2^)2 (2^)2 
and make there the following change of variables: 

(A.6) 
(A.7) 

R=^-^-'4 r2 = 2R + p2. {A.i 
Changing variables in Eq. (jA.Sp we obtain: 

4 f '^^P^ ^^^2 ^iR(A;i-gi-A:2)grpi(fci+fc2+2fc3-g2-g3)g»P2(g3-gi-g2) ^ 





rs - ri 




Pi 


2 


^3 




rs + ri 




P2 


2 


n 


r2 


ri rs 




2 


4 4 


r2 



(27r)2 (27r)2 (27r) 

Integration over and p2 gives three delta function and after the integration over gi, ki and (72 such 
that 

qi = q3 - k2 - ks (A. 10) 

ki = qs- ks (A.ll) 
92 = ^2 + fcs (A. 12) 

we obtain Eq. Q. 
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Appendix B: 

Let us consider together Eq. (fill) and Eq. (fTOl) : 

Si = J d^k ( Vi Vl ^Hy, k)) ^{y, k) ^{y, k) , (B.l) 

and 

The operator , because the rotational invariance of the problem (zero transfer momentum), can be 
rewritten in the following form: 

here we changed the variables: fc^ — > x. The whole expression Eq. (jB.ip now may be rewritten in the 
terms of variable x: 

Si = All f dx ^(x ^) Vl ^\x)) ^{x) ^{x) . (B.5) 
J ox ox 

Now we use , that 

f\x) oc x^Vl^\x) (B.6) 

(^(x) oc r dx l^\og{-). (B.7) 

Jx X ^ X 

Inserting these expressions in the Eq. (fB3]l we obtain: 

5, oc I dx^{x I-) (/t(x)) r dx ^ log(-) r dx" ^ log(-) . (B.8) 

J O X O X Jx X^ X Jx X ^ X 

Integrating Eq. (IB.SP by parts and omitting boundary terms of integration we obtain: 



5. a 2 / dx{^j\x)) [x f dx" ^ log(:^) dx' j . (B.9) 
The result of the second integration by parts is the following: 

5, oc -2 I dx/t(x) J^^ dx" log(^) f dx'^^+ (B.IO) 

+ 2 / dxfHx) (x r dx" r dx'^] + 



CO 



+ 2 / dx/t(x) ( X / - dx" log(^) l^ + fdx'i^^ 



that gives finally: 



Si oc 2 j —f\x) r dx' r dx' ^ + (B.ll) 

f dx , , /""^ , " fix) , ,x" , 
+ 2 / —f^{x)f{x) / dx Ars^log — . 

J X^ Jx X ^ X 

We see that this vertex is the same as in [8], see also ^6], excepting some not important for our consid- 
eration constant in front of the expression. 
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Appendix C: 



Let us consider expression Eq. 



+ (//-^(ri,r3)$(y,ri,r2)^>(y,r2,r3)) ( Li3$^(y, n, ^>^(y, r2, r3)) } . (C.l) 



and let us rewrite this expression in the foUowing form: 

/ f d?ri dPr'i f cfr'o d'^rl 



S.P = -\j dydyj / '-^^Hy,ny2)<^Hy,rl,r^ 



12' 2" 3 



IT 



P d^r-2 {(Li3«>(y,ri,r2)$(y,r2,r3))F-i(ri,r3) + 

^12 ^23 



+ (i? ^(ri,r3)$(y,ri,r2)$(y,r2,r3)) Li3}(5(y-y')(5^(^2 -^2) (C.2) 

From this expression for the Sj^p part of the action it is clear, that auto correlator Eq. (|3ip will be 
changed now to the following form: 



< i^[y,n,r2) ,ip{y ,r2,r3) >= (L12" ^>(2/,ri,r2))5(y - y )d (r2 - 
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2 ^2 



TT 



/"VIp ^2^^ |(Li3$(y,ri,r2)$(y,r2,r3))i/-i(ri,r3) + 

'^12 ^23 



+ (i/ Hn,r3)$(y,ri,r2)$(y,r2,r3)) Li3}5(y-yV^(^2 -^^2), (C.3) 
whereas the form of equation of motion Eq. (j29p will stay unchanged. 
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